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scaling with a length scale L(t) ~ t 1 / 2 . In particular, the autocorrelation function decays 
as We illustrate these properties by solving exactly the kinetic Potts model in 

d = 1. We then analyze a Langevin equation of an appropriate field theory to compute 
these correlation functions for general q and d. We establish a correspondence between 
the two-point correlations of the g-state Potts model and those of a kinetic Ising model 
evolving with a fixed magnetization (2/q — l). The dynamics of this Ising model is solved 
exactly in the large q limit, and in the limit of a large number of components n for the 
order parameter. For general q and in any dimension, we introduce a Gaussian closure 
approximation and calculate within this approximation the scaling functions and the 
exponent X(q). These are in good agreement with the direct numerical simulations of the 
Potts model as well as the kinetic Ising model with fixed magnetization. We also discuss 
the existing and possible experimental realizations of these models. 
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I. Introduction. 



Coarsening of domains of equilibrium ordered phases, following a quench from the 
disordered homogeneous phase to a regime where the system develops long range order, is 
widely observed in many physical systems such as binary alloys, liquid crystals, magnetic 
bubbles, Langmuir films and soap bubbles [1][2]. After the quench, domains of the 
different ordered phases form and grow with time as the system attains local equilibrium 
on larger and larger length scales. A dynamic scaling hypothesis suggests that at late 
times the system is left with a single length scale (linear size of a typical domain) which 
grows with time as L(t) ~ t n , where n depends on the conservation laws satisfied by 
the dynamics [1][2]. For systems with only two types of ordered phases (such as a 
binary alloy or the Ising model), the nonequilibrium coarsening dynamics have been 
extensively studied experimentally, numerically and by approximate analytical methods. 
Comparatively, much less is known when there are more than two ordered phases. A 
particular example of the latter class of models is the (/-state Potts model which has 
q ordered phases [3]. For q = 2, this corresponds to the Ising model, and there are 
experimental realizations also for q = 3,4, oo [3]. As q increases, the morphology of the 
coarsening patterns changes from one of large connected interpenetrating domains to one 
of more and more isotropic droplets. The limit q — > oo is known to correctly describe the 
evolution of a dry soap froth and the growth of metallic grains [4] . The reason for taking 
q = oo is to prevent bubbles represented by different Potts indices from coalescing. In a 
way, a finite but large q Potts model describes a froth in which wall breakage occurs with 
a probability of order 1/q (Fig. 1). Most studies of soap bubbles have so far focussed 
on geometrical properties of the froth. For instance, mean-field treatments [5] [6] and 
numerical simulations [4], as well as experiments [7], have addressed issues such as the 
joint distribution of bubble area and coordination number. However, while in other 
physical problems of coarsening the dynamic scaling and the correlation functions have 
been of central interest, no theoretical analysis exists of spatial and temporal correlations 
in a froth. 

2 



For the Potts model with finite q, there have been only numerical studies of the 
growth law of domains, confirming L(t) ~ t 1 / 2 [4] [8], and substantiating the scaling of 
the equal-time correlation function [8]. However, to our knowledge, the two-time corre- 
lation has never been studied previously, not even numerically. In this paper, following 
our recent letter [21], we show that the autocorrelation A(t) = (<j>(r, O)0(r, i)), measuring 
the correlation of the order parameter field 0(r, t) with its initial value 0(r, 0), decays 
as A(t) ~ L(t)~^( q \ Thus, while the growth exponent remains 1/2 for all q, the auto- 
correlation exponent depends explicitly on q. In [21], we established a correspondence 
between the dynamics of the g-state Potts model and that of a kinetic Ising model evolv- 
ing with a fixed magnetization (2/q — 1). In the present work, we explore more of this 
correspondence and derive several new results. 

For q = 2, the Potts model is by definition identical to the Ising model. For q > 2, 
our correspondence makes contact between the kinetic Potts model and a kinetic Ising 
model evolving with a fixed magnetization. This is useful not only because Ising models 
are conceptually much easier to think about, but also due to the fact that they are 
much easier to access experimentally, especially for the measurement of nonequilibrium 
dynamical properties. For example, in a recent experiment using video microscopy [27], 
the autocorrelation exponent has been measured for a liquid crystal sample confined 
between two parallel plates. The plates were designed such that molecules were bound 
to align along two perpendicular directions at the surface of the two plates. The two 
directions of the "helix" described by the molecules between the two plates (clockwise 
or anticlockwise) then represent the two possible "spin" states. This system then obeys 
Ising symmetry in d = 2 and corresponds to q = 2 (fixed 50%-50% mixture of "up" 
and "down" phases). The experimentally measured autocorrelation exponent \e xp = 
1.246 ±0.079 [27] is in good agreement with the numerical simulation of the d = 2 kinetic 
Ising model, X]y um = 1.25 ± 0.01 [9] (see also section VI). In Ref. [9], it was also argued 
heuristically that A = 5/4for<i = 2. The Gaussian closure approximation for q = 2 
leads to \qCA ~ 1-286 [15] (see also section V). In principle, for other values of q, X(q) 
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could be measured in such a system if one succeeds in maintaining a fixed concentration 
of "up" phase, different from 50%. 

For q > 2, exploiting our correspondence to an Ising model, another possible experi- 
mental system for measuring X(q) might be magnetic bubbles [26] [28]. Magnetic bubbles 
are a natural realization of the kinetic Ising model where increasing the magnetic field 
leads to the coarsening of the system. In principle, X(q) could be measured if a constant 
magnetization path in the phase diagram is chosen [28] . Note that in order to determine 
X(q), the "real time" scale which, as already mentioned, is mapped onto a function of 
the magnetic field in this system, is not explicitly needed since the definition of X(q) only 
involves the domains length scale, L{t). 

The paper is organized as follows. In section II, we solve exactly the dynamics of 
the d = 1 g-state Potts model and compute the equal-time and the two-time correlation 
functions. We find X(q, d = 1) = 1 for all q, and show that this problem is equivalent to 
the d = 1 Glauber dynamics for the Ising model with constant magnetization mg = 2/q — 
1. In section III, we present an exact theory of the coarsening of this Ising model with 
globally conserved magnetization mg = 2/q—l, in the q — > oo limit, and in any dimension 
d > 1. The distribution of droplet radii and the equal-time correlation functions are 
computed exactly and we find X(q = oo, d) = d. Section IV is devoted to the exact 
solution of an 0(n) model evolving with a globally conserved magnetization, in the 
n — > oo limit. The correlation functions are again calculated and X(q,n = oo, d) = d/2. 
In section V, we first establish a general correspondence between the dynamics of the 
g-state Potts model and that of the Ising model evolving with a fixed magnetization 
2/q — l, and then present a non trivial extension of the GCA [15] to this Ising model (a 
brief version of which was communicated earlier [21]). This approach gives very accurate 
results for general q and also reproduces the exact results obtained in the different limits 
mentioned above. Finally, in section VI, we present extensive numerical simulations for 
the Potts and Ising models separately, which confirm the mentioned correspondence and 
show that the GCA is indeed remarkably accurate. 
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As a prelude to the following sections, we first define the Potts model Hamiltonian 

as [3], 

« = -E^ ( L1 ) 

<ij> 

(jj takes integer values = l,...,q. The sum is over nearest neighbors, but can be 
extended to the next shell of neighbors, as is needed in T = Monte Carlo simulations, 
in order to avoid pinning to the lattice (see section VI). 

II. q-state Potts model in one dimension. 

We now consider the zero temperature Glauber dynamics of the g-state Potts model 
in one dimension. Contrary to higher dimensions where the domain growth is driven 
by interfacial tension, coarsening in d = 1 occurs via the diffusion and annihilation of 
kinks (Fig. 2-a). We start with a totally random initial condition where each of the q 
phases is present with equal density c = \jq. The zero temperature dynamics proceeds 
as follows: a spin is selected at random and its value is changed to that of either of its 
neighbors with equal probability. It gives rise to three possible situations as shown in 
Fig. 2. We now focus on one particular phase with an arbitrary Potts index a = I and 
define an indicator function <f>(x, t) which is 1 if the site x is occupied by this l-th phase 
and otherwise. In this new two-phase system, the density of l's is c = 1/q, and the 
O's, representing the (q — 1) other phases, have a density 1 — c. For later convenience, 
we prefer the <f> = 0, 1 convention to the more usual spin representation <p = ±1- The 
density c is related to the usual magnetization mg by the relation tuq = 2c — 1. Then 
the dynamics of <j) is governed by the Glauber dynamics [10] of the Ising model with a 
constant magnetization uiq = 2/q — 1. This is illustrated in Fig. 2-b. 

Using this mapping, the equal-time correlation function G(r,t) = (<j>(r, t)<f)(0, £)), in 
the scaling limit, is given by [10], 

5GM)= (l-I)erfc(^) + i (III) 

where erfc(x) = J^?°exp(— v?) du. We now point out that in the q — > oo limit this 
correlation function is related to the probability distribution of spacings between do- 
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main walls (kinks in d = 1). The quantity qG(r,t) measures the probability that the 
spins at and r have the same value at time t. For finite q, this is different from the 
probability that they belong to the same domain. However, when q — > oo, two spins 
which are equal are necessarily in the same domain (bubble). Then, qG(r,t) is the 
probability of having a domain of length r or more. It is then a standard result that 
the normalized spacing distribution of kinks P(r,t) is given by P(r,t) = 2qL(t) — ^ ' . 
For large times, this distribution also obeys scaling: P(r,t) = L(t) p(r/L(t)), where 



obtained by Derrida et al. by using an analogy to a random walk problem [11]. Interest- 
ingly, this Wigner distribution is also identical to the eigenvalue spacing distribution for 
a real random matrix [12]. 

We also note that this result coincides with the spacing distribution in the reaction 
diffusion model A + A — > A [13]. This can be understood in the following way. Repre- 
senting a kink by a particle A on a d = 1 line, the dynamics of the Potts model can be 
mapped to that of a reaction diffusion model where the particles A diffuse, annihilate 
and coagulate according to the following rules: each particle undergoes diffusion until 
two of them meet in which case they either annihilate (A + A — > 0) with probability 
l/(q — 1), or coagulate (A + A — > A) with probability (q — 2)/(q — 1) (see Fig. 2-c). For 
5 = 2, the particles only annihilate and hence, this is equivalent to the Glauber model 
[10], whereas in the q — > oo limit they only coagulate. Both these limits have been stud- 
ied previously [13] [14]. In the Glauber case, the two-point correlation function has been 
calculated analytically but the spacing distribution is still unknown. On the other hand, 
for the A + A — > A model, the spacing distribution is known exactly but there was no 
analogy to any spin model and hence correlation function. Our present result establishes 
that the dynamics of these two problems are two different solvable limits of that of the 
g-state Potts model. 

We now consider the two-time correlation function C(r,t) = (6 a r r j-\ ,7(0,0) ~~ V? 2 )- 
In terms of the indicator field </>, C(r,t) = (<p(r,t)<p(0,0) — 1/g 2 ), which is expected 
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This result, in the q 



+00 limit, was first 



to scale as C(r,t) ~ L(t)~ x c(r / L(t)) [9]. Thus, the autocorrelation A(t) = C(0,t) 
decays as A(t) ~ L(t)~ x , where A is a non trivial non equilibrium exponent [9]. In 
d = 1 and for all q, we find that C(r, t) satisfies a diffusion equation and is given by 
C(r,£) ~ t _1 / 2 exp(— r 2 /2t). Since the length scale L(t) ~ t 1 / 2 , we establish the scaling 
of C(r,t) and find A = 1 for all q. This is consistent with our general exact result that 
A = d in the q — > oo limit, as we argue in the following sections. 

In one dimension, we have shown that the evolution of the two-point correlation 
functions in the kinetic Potts model can be exactly mapped to that of the Ising model 
with fixed magnetization (2/q — l). In section V, we will show that this correspondence 
essentially holds even in higher dimensions. This fact motivates us to study the dynamics 
of the Ising model evolving with fixed magnetization. In the following section, we show 
that this problem can be exactly solved in the q — > oo limit. 

III. Small volume fraction limit. 

In this section, we study the coarsening of a magnetic system with globally conserved 
magnetization mg = (2/q — 1), or equivalently with a density c = \jq of the minority 
phase (up), in the c — > limit. This becomes a modified version (suited for globally fixed 
magnetization) of the celebrated Lifshitz-Slyozov (LS) theory [16] [25] which describes the 
coarsening of a two-phase system with local conservation (model B [17]). On a discrete 
lattice, both models are described by spin-exchange Kawasaki dynamics. However, in one 
case the exchange occurs between spins on any two arbitrary sites (global conservation) 
whereas in the other case, the exchange occurs only between two sites which are nearest 
neighbor to each other (local conservation). For c = 1/2, the globally conserved model 
(hereafter called model A(c)) has been shown to be in the same universality class as the 
standard model A dynamics, so far as the dynamical exponent and the domain growth 
exponent are concerned [18] [19]. In section VI, we will show that the autocorrelation 
exponent A is also the same for a strict global c = 1/2 conservation and for model A, 
namely, A ~ 1.25. However, for c < 1/2, we find (see section VI) that A is explicitly a 
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function of c. At the end of this section, we show (within the framework of LS theory 
which is exact in the zero volume fraction limit) that lim c ^o A(c) = d. We recover this 
exact result also from our approximate treatment of the appropriate field theory (see 
section V). This result is further verified from the direct numerical simulations (section 
VI) of the g-state Potts model (large q limit) as well as the model A(c). 

In the c — > limit, the coarsening pattern consists of circular bubbles which are 
growing, but are always far apart from each other so that they never coalesce. A mean- 
field treatment assuming no correlation between these bubbles should then be exact in 
the limit c — > 0. In the following, we determine N(R,t), the density of bubbles with 
radius r at time t. 

The equation of motion for the radius R{(t) of bubble i is, 

% = -lQ + Xi (IIU) 

The first term of the right hand side, in which the surface tension dependent coefficient 
has been normalized to 1, contributes to make the bubble shrink, in order to minimize 
locally the interface length between the two phases. The second term, Ai(t), is an effective 
time-dependent magnetic field, playing the role of a Lagrange multiplier (see section V) 
fixing the constraint that J2i Rf must be a constant proportional to cV, where V is the 
volume. Note that this equation of motion Eq. (III.l) can also be directly derived from 
the model A equation of motion with this additional Lagrange multiplier needed to keep 
the magnetization fixed. The density N(R, t) satisfies a continuity equation of the form, 

dt dR\dt J y 1 

At late times we look for a scaling solution of the form N(R,t) « L(t)-( d+ ^F(R/L(t)), 
where L(t) is proportional to the average radius of a growing bubble. The conservation 
law demands that the prefactor decays as L(t)~( d+1 \ Inserting Eq. (III.l) and the 
scaling form in Eq. (III. 2), one easily sees that all terms are of the same order provided 
is a constant, so that L(t) ~ t 1 / 2 . Thus, in contrast to model B dynamics where 
L(t) ~ t 1 / 3 , the model A(c) has the same growth law as model A. We now set = 1 



and X\(t) = v/L(t), where v is to be determined by imposing that the solution of the 

scaling equation has physical limits. This form for Ai is justified by the fact that, in 

Eq. (III. 1) , Ai scales as R . This fact will be physically justified in section V in terms 

of balance between the interface and magnetic energies of a droplet. From now on, we 

present explicit results for d = 2, but the generalization to d > 2 is straightforward. The 

scaling function F then satisfies the following differential equation: 

1 _ ^ r 2 

F\x) = , 2 - F(x) (III.3) 

x(x 2 -ux + 1) 

The condition that F(x) goes to zero for large x demands that v = 2, so that the 
function multiplying F(x) in Eq. (III. 3) has a double pole (at x = 1). Then, solving this 
differential equation we get, 

x ( 2 \ 

Fix) = A- -j exp ( | , for x < 1 

1 ' (1-x) 4 I 1-J " (III.4) 



= 0, for x > 1 

where the constant A is determined from the conservation condition. This scaling func- 
tion is plotted in Fig. 3 and appears to be wider than the LS form for model B (local 
conservation) [24]. 

We now compute, in the limit c — > 0, the equal-time correlation function G(r, t) = 
{4>(r : t)(f)(0, t))/c, where <fi is the density field as defined in the previous section. The func- 
tion G(r, t) can be computed from the normalized radius distribution function Nq(R, t) = 
N(R, t)/J iru 2 N(u, t) du = N(R, t) /cV in the following way. By definition, 

^ v ' i,3 

where the indices i and j run over all the bubbles and Xi( x ) is the characteristic function 
of the z-th bubble. In the c — > limit, bubbles are strictly circular and Xi( x ) = &(Ri — l x l) ; 
where 9 is the usual step function and Rj is the radius of the i-th bubble. In this low 
area fraction limit, the bubbles are far apart from each other and for finite r, only the 
terms corresponding to i = j contribute to the sum in Eq. (III. 5). Then, we obtain the 
exact result in the c — > limit: 

G(r, t) = j dRNo(R, t) J rf 2 x Xi? (x + r) Xii (x) (III.6) 



The second integral is just the overlap area between two disks with radius R, with their 
centers separated by a distance r. The final expression is, 



G(M) = 2 j^dHJftVbW*) (arccos (JL) - JLjl - J^j (III.7) 

In the scaling limit, No(R,t) « ^^r S F(R/L(t)), where is given by Eq. (III.4), 

and therefore, from the above equation, we find explicitly that G(r,t) ~ g(R/L(t)). 
Thus, dynamic scaling is established. The function g(x) with x normalized such that 
g(l) = 1/2 is shown in Fig. 4, and is seen to be in nice agreement with large q Potts 
model simulations, and the field theory results presented below. Notice however that the 
mean-field equal-time correlation function has a finite support. 

In the c — > limit, it is also simple to calculate the autocorrelation exponent A. Since 
bubbles do not coalesce in this limit and their centers do not diffuse, the autocorrelation 
A(t) ~ (0(x, £)0(x, 0)) is essentially the survival probability of a bubble up to time t. 
After time t, the number of bubbles left is N(t) = J dRN(R, t) ~ L{t)~ d , and therefore 
A = d, in d dimensions. This result is confirmed by numerical simulations presented in 
section VI. 

IV. Large n calculation. 

Another example where the correlation functions (both equal-time and two-time) and 
the exponent A can be calculated in presence of a time dependent magnetic field Xi(t) (to 
keep the average magnetization uiq fixed), is the large n (n — > oo) limit of the classical 
0(n) vector model. The 0(n) model is described by an n component order parameter 
field 4> ( r 5 1) = [4>\ (r, t), . . . , </> n (r, £)] and a coarse-grained Landau-Ginzburg free energy 
functional, 



n(7) = lfd d r 



Jl ft, / — i.2 



2 



a 



(IV.l) 



where \\{t) is a Lagrange multiplier to keep y J2r 0a ( r )^) fixed at mo = 2/q — 1 where 
V is the system size. The model A equation describing the overdamped relaxation is in 
general nonlinear and hard to solve. However, in the large n limit, this equation can be 
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linearized in a self-consistent way [20], and therefore, is exactly solvable in that limit. In 
Fourier space, this linearized equation (at T = 0) reads, 

^jM = _[fc2 + rQ + uSo (t)]<p(K t) + Xi(t)W5 k=Q (IV.2) 

where self-consistency for the structure factor demands that, So(t) = y Y^h <5(k, £), where 
<S(k, i) = (0(k, t)<f>(— k, t)). Note that ro < 0, since we are in the ordered phase. We have 
dropped the subscript a from <f> a since in the large n limit, the different components 
of <fi are completely uncorrelated with each other. The Lagrange multiplier X\(t) is 
determined from Eq. (IV.2), by demanding that 0(k = 0, t) = moVV and is given by, 
X\(t) = mQ[ro+uSo(t)]. We use random initial conditions, for which S(k, 0) = A — ijiq for 
k ^ 0, where A is of 0(1) and S(k = 0, 0) = m%V . Writing, Q(t) = J *[r + uS (t')] dt', 
we get from Eq. (IV.2), for k ^ 0, 

5(k, t) = (A - m§) exp[-2(A; 2 t + Q(t))]. (IV.3) 



The self-consistency condition now reads, 

dQ 

1F = r + u^ , y 



^ = r + W m§ + ^E^(k,t) (IV-4) 



k 

Plugging in the expression for S(k,t) from Eq. (IV.3) into Eq. (IV. 4), and taking the 
thermodynamic limit V — > oo, we get 

^ = r + um\ + u(A - m§)H- rf / 2 exp[-2Q(t)] (IV.5) 

where T is a constant which depends on the dimension d and the upper cutoff A of the 
theory. This equation is solved by making the ansatz, Q(t) = A + B logt. The left hand 
side of Eq. (IV.5) then decays as £ _1 , whereas the leading order term on the right hand 
side is a constant. So, for consistency, one needs B = —d/4, so that the leading order 
term on the right hand side is identically zero. As a consequence, the structure factor 
for k ^ (Fourier transform of the equal-time correlation function) can be written in 
the form, 

S(k, t) w L(t)- d s(kL(t)), with 

(IV.6) 

L(t) ss t 1/2 , and s(x) = Cexp(-2x 2 ) 
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where C is a constant. We thus obtain the scaling of the correlation function with the 
expected domains length scale, L(t) ~ t 1 / 2 . 

Similarly, the two-time correlation function, C(k, t) = (</>(k, 0)(p(— k, t)), for k ^ 0, 
evolves as C(k, t) = (A — mjj) exp[— (k 2 t+Q(t))]. The autocorrelation function, defined as 
A(t) = (0(r, O)0(r, £)) - mg, then decays as t~ d ^ ~ L(t)~ d l 2 . Thus, the autocorrelation 
exponent is A = cZ/2, as in the mo = case [20]. This, however, is not unexpected because 
the limit n — > oo decouples the different k modes, and the time dependent magnetic field 
just couples to k = mode which has no effect on the evolution of the k ^ modes 
apart from modifying the prefactor. Therefore, in order to see the dependence of A on 
mo, one has to include the 0(l/n) corrections which is a very hard task [2]. 

V. Field theory and Gaussian closure approximation. 

In our earlier letter [21], we constructed a field theory of the g-state Potts model in 
terms of the coarse-grained "occupation density" fields {<^(r, t); I = 1, 2, . . . q} such that 
4>l assumes the value 1 in the interior of the l-th ordered phase and decays continuously to 

outside. Consequently, inside any "bubble" of one phase, only one of the 0/'s is close to 

1 and the others are all close to 0. We thus require a potential with q degenerate minima 
at [1, 0, 0, ... 0], [0, 1, 0, ... 0], ... [0, 0, 0, ... 1], which prevents two different bubbles from 
sharing the same position in space. A suitable free energy functional is [21], 

F[{4>,}\ = / d d r £, (i(V<« 2 + V(<h)) - Ai t f>, - lj 

+A2 (sH) : 



2 

q-1 



(v.i; 



where Ai(r, t) is a Lagrange multiplier enforcing the constraint J2i<t>i = 1) an d A2 is a 
constant of 0(1) such that the state [1/q, . . . ,1/q] is unstable, and V{4>) ~ </> 2 (l — 4>) 2 is 
the usual double well potential with minima at and 1. Then, the equation of motion 
is, 

f)rh, „ . / 1 \ T JL / ^\ 2 n -^ 

(V.2) 



001 V 2 o, V f (o,) ■ A 1 1A, L- 1 



dt \ <i 
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i'=l 



and Ai = | Ez V(<f>i) by demanding £ z <t>l = 1 in Eq. (V.2). Note that for g = 2 (Ising 
model), Ai = by virtue of the condition <f>i + (fi 2 = 1 and A2 can be chosen to be since 
this term only renormalizes V. Then, one recovers the usual time-dependent Landau 
equation for the Ising model. For q > 2, we also note that this evolution equation has 
a form similar to that of Eq. (2.10) of Lau et al. [8]. In Fig. 1, we show a late time 
configuration of domains generated by Eq. (V.2) for q = 30. 

The two-point correlation function for the g-state Potts model is defined as (2(12) = 
Y,f =1 (4>i(ri, £i)0/(r 2 , t 2 )) and therefore equals g(0/(r 1; £i)0/(r 2 , t 2 )) due to the symmetry 
between the q phases. Here, "12" is a short-hand notation for the pair of space-time 
points (ri, t\) and (r2, t 2 ). Due to the isotropy and translational invariance in space, the 
only spatial dependence of these correlation functions is through r = |ri — r 2 |. Denoting 
the equal-time correlation function [t\ = t 2 =t) by G(r,t), we get from Eq. (V.2): 
1 FtC 

Y^ = v2G ~ q t] ( y/( ^ (r ' t)] ~ Al(r ' 

- 4A 2 g ^(0, t) ^(r, t) - (j2 t) - 1 

Note that the two-time correlation function satisfies a similar equation. 

Our first approximation is to replace 4>i'^ by its average q{4>i^) = (2(0, t) in the 
third term on the right hand side of Eq. (V.3), which becomes exact in the q — > 00 
limit. Furthermore, the scaling solution G(r,t) = g(r/L(t)) must satisfy g(0) = 1, so 
that we can drop the term 4A 2 ((2(r, t) — l/q)(G(0,t) — 1) so produced. Thus, the third 
term, although important in the evolution of 0/ since it provides stability to the bubbles, 
is not crucial in the evolution of the correlation functions, at least in the scaling limit 
of the large q model, but also for q = 2 for which this term is simply absent. The two 
boundary conditions for G(r, t) are: (i) As r — > 0, G(r, t) — > 1 and (ii) As r — > 00, 
G(r,t)^q{(/>i(0,t)){</> l (T,t)) = l/q. 

Next, using E0/(r, t) = 1, we get (E <M r ) £)Ai(r, t)) = (Ai)(t), and then given the 
symmetry between the q phases, we can write (0/Ai) = (l/g)(Ai) = (0/)(Ai). Thus, 
without approximation, we replace in Eq. (V.3) Ai(r, t) by its average (Ai(r, t)), which 
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)> 



is simply a function of time. As a result, Eq. (V.3) reduces to an equation involving only 
a single field </>/(r, t): 

1 r)C 

Yq[ = V 2 G - q (0 Z (O, t) (v'iU*, t)) - (Ai))) (V.4) 
Interestingly, Eq. (V.4) is also the evolution equation for the two-point correlation in 
an Ising model evolving with fixed average magnetization (mo) = 2/q — 1 (equivalently, 
with a density of minority "up" spins fixed at 1/q). The droplet domains of the minority 
phase in this Ising model would correspond to the bubbles of a particular phase in the 
Potts model, with the majority phase corresponding to the remaining (q — 1) phases. 
(Ai) acts as a time-dependent magnetic field which prevents the minority phase from 
disappearing at T = 0, and keeps the magnetization constant. 

Thus, from now on, instead of the original Potts model, we consider the coarsening 
in the Ising model with globally conserved magnetization through a time-dependent 
magnetic field (Lagrange multiplier). The magnetization mo is related to the value of q 
of the Potts model via ttlq = 2/q — 1. 

The calculation of the two-point correlation function of this problem can be performed 
approximately by extending the Gaussian closure scheme as developed by Mazenko [15] 
for the usual Ising model where the magnetization remains fixed at mg = 0. The essence 
of this approximation scheme is to invoke an auxiliary field m(r, t) which is related to 
the order parameter field (f>(r,t) via a non linear transformation <f>(r,t) = a(m(r,t)). 
The idea is to find a field m(r, t) which varies smoothly across the interfaces as opposed 
to the original field 0(r, t) which changes abruptly from nearly to nearly 1 across an 
interface. So, the non trivial part of the scheme is to choose the appropriate mapping 
function a. In the case of the Ising model with zero magnetization, Mazenko argued that 
the function a should be chosen as the equilibrium profile of <fi near an interface, which 
is determined by the solution of 

S = ^ (v ' B » 

with the boundary conditions a(m) — > 1 as m — > +oo and a(m) — > as m — > — oo. 
The solution is a(m) = [1 + tanh(/?m)]/2, which at late times can be replaced by a step 
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function, since domains grow (with L(t) ~ t 1 / 2 ) whereas the interface width (~ 
related to the coupling constant in front of V) remains bounded. The auxiliary field 
m(r, t) can be interpreted as the distance from the nearest interface. The next part of 
the approximation is to assume that m(r, t), being smooth across the interface, has a 
Gaussian distribution. The virtue of this "minimal" approximation is that it facilitates 
analytical calculation yielding non trivial results for the correlation functions and the 
exponents which are in good agreement with simulations, at least in the nonconserved 
case. 

In order to extend this approximation to our problem with globally conserved mag- 
netization, we start off with the assumption that there exists, as in the nonconserved 
Ising nonlinear transformation (j)(r,t) = a(m(r,t),t) with a Gaussian auxil- 

iary field. However, several important modifications need to be done in carrying out 
this extension from the simple Ising case. First, (0(r, £)) must be strictly fixed at 1/q 
at all times (by virtue of the time-dependent field (Ai)(t)), as opposed to the Ising 
case where, for a critical quench, ((f)) = 1/2 automatically. This also necessarily im- 
plies that the mean of the distribution of m is nonzero. The first and second mo- 
ments of the Gaussian distribution, (m(r, £)) = fh(t), ({m(r,t) — m(t)} 2 ) = Co(t) are 
space independent due to translational invariance. The complete correlation function 
({m(ri,ti) — ra(£i)}{m(r2,t2) — m (^2)}) = C(12) must be determined self-consistently 
as is C(12) in the Ising case [15]. Now, from the condition (0(r, £)) = 1/q, we get 
^2^(7 ^-^oo <J { m ) exp[— (m — m) 2 /2Co] dm = 1/q. Replacing a(m, t) at late times by the 
step function 6{m) (which is 1 for m > and for m < 0), and thereby neglecting 
terms that are of lower order in t, we get m(t) = — \J 2Cq (t)erfc~ 1 (2/q) . Note that for 
q = 2, rh = as expected. For later convenience, let us also define the correlation func- 



tion /(12) = C(12)/ v / Co(l)C (2) and denote it by f(r,t) when t\ = t 2 = t. Note that 
/(0, t) = 1, and / — > as r — > oo. Even for t\ = we well keep on writing Cq(1) and 
Co (2) explicitly, although this two numbers are equal, since we will use formal derivatives 
with respect to Cq(1). 
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The second important difference from the simple Ising case is the choice of the map- 
ping function a(m(r, t), t). The explicit time dependence introduced via (Ai) modi- 
fies the local equilibrium profile, thereby precluding the choice of the stationary profile 
[1 +tanh(/?m)]/2. In fact, since the mean of the field m(r, t) is time-dependent (through 
Co(t) which scales linearly with t as we argue below), one can expect to get a "sig- 
moid" shaped solution only in a moving frame, with a velocity suitably determined to 
neutralize the time-dependence introduced via (Ai). Thus, making the transformation, 
r' = r + a(t)h where n is an arbitrary unit vector and a(t) is to be determined, and 
demanding an equilibrium solution, i.e., ^ = to leading order in time, we find the 
appropriate equation for a(m,t): 

do da do tW/ . . , TJ . „. 

We now fix a(t) from the condition that the average value on both sides of Eq. (V.6) 
should be identical. The average on the right hand side is zero by definition of (Ai). 
Then from Eq. (V.6), we get, 

da_ (o"(m)) _ gg du exp ~^fo"(u) 
dt (^H> f+£duexp-k^jfa'(u) 

^ ow ^^77?' * / ($m) ls ca l cu l a ted replacing a(m, t) by 9(m) at late times, so that a'(m) ~ 
8{m) and a"(m) xs 8'{m). This yields, 

da d , / m 2 \ rh _ 

^"-^ log ( exp "^J = Q (v - 8) 

Anticipating a scaling solution, we find that a(t) ~ \JCo(t) ~ L(t), which is expected, 
since L(t) is physically the only remaining length scale at late times. From Eq. (V.8) 
and the expression of m as a function of Co(t), we see that (Ai) ~ L(t) -1 , which can be 
understood on physical grounds: local equilibrium of a bubble and its interface requires 
that the surface tension energy Eg ~ L(t) d ~ l should balance the magnetic energy Em ~ 
{\l)L(t) d (see also section III). 
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Using the fact that m(r, t) has a Gaussian distribution, the correlation function 
G(r, t) is given by, 



Q1 f 

G(r,t) = , = / dx\ dx2 <j{x\)a{x2) x 

27r,/C (l)C (2) ■/ 



exp 



7 2 / (>i - mi) 2 (x 2 - m 2 ) 2 , _ v _ . / 

„ 7T\ 1 7777^ h2 ai-mi x 2 -m 2 , -— 

Co(l) Co(2) ^/Co(l)C (2) # 



(V.9) 



where 7 = 1/^1 - f 2 and we recall that /(12) = C(12)/^C (1)C (2) and that the 
arguments "1" and "2" denote (r, i) and (0, t). The derivative with respect to m involved 
when inserting Eq. (V.6) in Eq. (V.4) are more easily expressed in the Fourier space 
associated to the variable x. G(r, t) then takes the form, 

q r 

G(r,t) = J dkidk2<j(ki)a(k2)x 
k 2 ^ k 2 ^ ( V - 10 ) 



exp 



^Cq(I) - -fC (2) - k x fc 2 C(12) + %mi + ik 2 m 2 



where a is the Fourier transform of a. From this expression one finds that, 

/ ^rfV^A n dG 
^( 2 )X72( 1 )) = 2: 



dm 2K 7 aC (l) 
/ ^dtr^A dG 



mi 



(v.n) 



dm / dmi |c (i) 



Using this result and noting that for large time q^r\) = 2 g^|^ , the second term on the 
right hand side of Eq. (V.4) can be written as, 

\ (V.12) 
<9G ami dG )_ 2 9G f dG 



dC (l) lrfll dC {l)dfh im i)J dC (l) C (l)df 
Therefore, the evolution equation for the correlation function can be expressed as, 

5§H 2G + a7)W> (v ' 13) 

where Q(f) = /^r- Interestingly, we notice that the last equation is identical in form 
to that of the Mazenko equation [15] [21] for the critical Ising case with the exception 
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that G(f) has different expressions in the two cases. However, this seems accidental 
because in our problem we need to invoke a moving frame and therefore a different 
profile function satisfying Eq. (V.6). A naive application of the Mazenko theory with the 
mapping function determined by the equilibrium profile Eq. (V.5) leads to inconsistent 
and unphysical results. For q = 2 (c = 1/2, Ising critical), the velocity of the moving 
frame ~j/jr ~ "VCq is zero identically, and our expression then reduces to the Mazenko 



Replacing a by the 9 function in Eq. (V.9), we get the leading term for G(f) for 
large time: 



where p = erfc (2/g). In principle, the function G(f) can be inverted, so that Q 
is implicitly a function of G. Note that this function has the correct short (as / — > 1, 
G(f) — > 1) and long distance (as / — > 0, G(f) — > 1/q) behaviors. In addition, for q = 2, it 



reduces to the Ising case [15], G(f) = ^ tan [y (1 + /)/(l — /)]. We now substitute the 
scaling form G{r,t) = g(r/L(t)) in Eq. (V.13). A scaling solution is obtained provided 
Co(t) ~ L(t) 2 ~ t, which leads to the expected form for L(t). The fact that Co(t) is 
proportional to L^(t) is consistent with the definition of Co(t) as a two-point correlation 
function of the field m, which has the physical meaning of a distance. The condition 
L 2 (t) ~ t is obtained by demanding that all terms in Eq. (V.13) obtained by plugging 
in a scaling form for the correlation function are of the same order as a function of time. 
More precisely, setting Co(t) ~ At/ \i and x = r/L(t), we get from Eq. (V.13), 



which defines a closed eigenvalue equation for the scaling function g. The eigenvalue 
has to be determined by matching the short and long distance behaviors of g(x). The 
autocorrelation exponent A is then related to fx via the relation A = d — fi/2, following 
an argument due to Mazenko [15] that we adapt to our problem below. 



result. 




(V.14) 





(V.15) 
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Following the same line of arguments as used to derive the evolution equation for G, 
we find that the two-time correlation function C(r,t) = q((p(r,t)<p(0,0)) — q^ 1 satisfies 
the equation, 

where Q(f) = /§^, and C(f) has the same / dependence as G(f). Since the two- 
time correlation function decays with time, its value is very small at late times and 
then Q ~ C. With Co(t) ~ At/n, this linear equation can be solved, and C(r, t) ~ 
t -(d/-n/2)/2 exp (_ r 2 / 2 £). Therefore, the autocorrelation A(t) = C(0,t) ~ H d -M/2)/2 _ 
L ( t )-(d-/x/2) j and x = d-/i/2. 

In the Q — > oo limit, it is possible to solve Eq. (V.15) analytically. Neglecting terms of 
0{l/p^) and using erfc(p) ~ exp(— p^)/p\^~ for large p, we find g(f) ~ ^^erfc 



Note that as / - 0, g(f) = 1/q + /(l + 2p 2 )/ ? + 0(/ 2 ), and as / - 1, g(f) - 
1 — P\j2(l — f)/n. Then, expressing the function Q in terms of g itself, we find essentially 
three regimes. As g — > 1/q (large distance), Q(g) ~ g — 1/q, and as g — > 1 (short 

2 

distance), Q(<7) ~ ^fz^y- I* 1 the intermediate regime, <?* <C g <C 1 where g* ~ log(q)/q, 
Q(g) ps p 2 g/2. Note that, as g becomes larger and larger, the window over which Q((7) 
behaves as g — 1/q becomes smaller and smaller and for a large range of distance one has 
/ ps 1. First consider the small x behavior of Eq. (V.15). Using Q(g) ~ p 2 /ir(l — g), 
we find that g(x) — > 1 — P^J n ^_-^ x ^ reflecting the presence of sharp interfaces. This 
reproduces Porod's law [22], namely that the structure factor scaling function F(y) : the 
Fourier transform of g(x) : decays as for large argument y. 

For very large g, since in the interesting range of distance one has /(r, t) ~ 1, one can 
find a very simple parametrization for g and Q(g). Using the new variable u ps P\f^r~-> 
we obtain, 

g{u) = erfc(w) 

n< \ 2 p^ P 2 2 (V - 17) 

Q{u) =p K(u) = !— exp — u 
2^/iru 

Because of the overall factor p 1 in the expression of Q(u), and since p grows with q 
(p 2 ps log((/)), we expect v = lim^oo n{<l)p* to be finite. Then, using the parametrization 
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of Eq. (V.17), u{x) satisfies the following eigenequation: 

„" _ 2uu' 2 + u' (— + x )-^ = (V.18) 

Porod's law gives u ~ ^ 2(d-l) x ^ or x ~^ ^' anc ^ ^ e ^ ar S e x limit is easily found to be 
it x/y/2. By matching the two regimes, one can check that it (a;) = x/\/2 is a solution 
of Eq. (V.15), provided the eigenvalue v satisfies v = 2(d — 1). The eigenfunction is 
then g(x) = erfc(x/v / 2). Moreover, for large q, \i ~ 2(d — l)/p 2 , where ~ log((/), and 
therefore A ps d — (d — l)/logg. 

We note that, for d = 1, the eigenvalue problem can be solved directly for any q. The 
small x behavior of g(x) obtained from Eq. (V.15) implies \i = 0. The scaling function 
then satisfies the differential equation g" + xg' = 0. As a consequence, g{x) coincides 
with the exact solution (Eq. (II. 1)) of the d = 1 (/-state kinetic Potts model presented in 
section II. \i = leads to A = 1, which was also a result of the exact d = 1 calculation. 
Thus, the Gaussian closure approximation is exact in d = 1. 

VI. Numerical simulations. 

We now compare our results with the direct T = simulation of the (/-state Potts 
model. We have also simulated directly our field theory (Eq. (V.2)) and found that 
it evolves in a similar way as the Potts model (see Fig. 1) with domains growing as 
L 2 (t) ~ t. However, in the field theory, one needs as many fields (4 bits real) as Potts 
indices which limits the maximum lattice size (~ 120 x 120), q (q ma x ~ 50) and the 
number of time steps. The determination of A requires large lattice (especially for large 
q) and large numbers of MC steps (typically 10 4 or more) which is easier to achieve 
in the Potts model simulation. The calculations have been carried out at T = on 
a 800 x 800 square lattice with equal coupling to nearest and next nearest neighbors 
(NNN). NNN interactions are needed at T = to avoid pinning to the lattice for q > 2. 
It also ensures a better isotropy of surface tension and thus of bubbles. We optimized 
the MC procedure by only selecting surface Potts spins for updating, since only these 
can be flipped at T = 0. For q < 20, the results for 30 to 40 samples were averaged 
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(more than in [21]), whereas 10 to 20 samples were found to be sufficient for q > 20, due 
to smaller fluctuations for L(t) and A{t) with increasing q, yielding, to our knowledge, 
the most extensive simulations to date. 

We found L 2 (t) ~ t for all q, confirming a result already obtained in previous studies 
[4] [8]. We also observed good scaling of the correlation functions. In Fig. 4, we compare 
the equal-time correlation function for q — > oo as given by the Potts model simulation, the 
direct simulation of the field theory for q = 50 fields on a 120 x 120 lattice, the mean-field 
theory of section III, and by our approximate theory, and find good agreement between 
the numerical results and the two theoretical ones. For soap bubbles, g(x) measures the 
probability that the point x belongs to the same bubble as the origin. Since bubbles 
are essentially isotropic, we expect g to be closely related to the distribution of radii. 
Since we find that g(x) has a Gaussian tail, the distribution of areas (A ~ x 2 ) has a 
Poissonian tail. This result is consistent with a maximum entropy [23] and a mean-field 
[5] [6] calculation. More precisely, the second derivative of g{x) is proportional to the 
interface spacing distribution on a linear cut through the froth, which has a Wigner form 
[12] in our case. 

A more interesting test of our theory concerns the computation of the autocorrelation 
exponent A. The mean-field result of section III predicts A = d, whereas the large n 
calculation (section IV) leads to A = d/2. The more sophisticated Gaussian closure 
approximation of section V gives less trivial and more accurate results: in Tab. 1, 
we present the values of A generated in Potts simulations, and compare them to those 
obtained from the (numerical) solution of the one-dimensional eigenvalue problem of Eq. 
(V.15). We find a reasonably good agreement. Notice that for large q (q > 100) the 
exponent A obtained from Monte Carlo simulations is probably slightly overestimated. 
This is due to the fact that for a finite lattice q = oo is actually realized for a finite value of 
q, so that one can expect that the effective number of different phases is q e f f > q. Finite 
size scaling indeed confirms this fact, although X(N) does not seem to have a simple 
form. In fact, for the values of q presented here, the finite size correction is comparable 
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to the statistical error bars. For instance, for q = 200 and a 300 x 300 lattice, we find 
A = 1.84 ± 0.02 instead of A = 1.82 ± 0.02 for a 800 x 800 lattice. We also note that 
A from the simulation saturates very slowly to its g — > oo value, as predicted by our 
asymptotic result. For soap bubbles, q = oo and A = 2 (A = d in dimension d). Indeed, 
the choice q = oo eliminates the coalescence of bubbles with identical index, so that 
A(t) ~ N(t)~^ ~ L(t)~ d , where N(t) is the number of remaining bubbles at time t (see 
also section III) . 

We also tested the correspondence between the two-point correlation functions of 
the g-state Potts model and that of the Ising model with a fixed average magnetization 
(mo) = 2/q — 1, as suggested by our analysis of section V. The dynamics is described by 
an infinite range Kawasaki dynamics where two randomly selected spins are exchanged 
with probability 1 if the energy is lowered, and with probability 1/2 if the energy is not 
changed [18]. Again, the algorithm is optimized by keeping track of the movable (surface) 
spins and by rescaling properly the unit of time. This last aspect is not important for 
the determination of A which does not involve the real time explicitly, but only L(t). 
For q > 2, this simulation appears to be much more delicate than for the g-state Potts 
model or for the Ising model with q = 2. Indeed, the choice of the initial conditions is 
crucial as it is already known [24] in simulations of the dynamics following an off-critical 
quench in the conserved order parameter (model B) case. For instance, one can start 
from an assembly of preformed circular bubbles with a distribution of radius given by 
the mean-field expression of Eq. (III. 4) . The bubbles are then randomly placed on the 
lattice. Although attractive, this procedure leads to very long transient times since the 
correct correlations between bubbles are long to establish through merging of bubbles. 
In other words, the (absence of) correlations introduced in the initial state are very long 
to destroy, a phenomenon which is amplified for large q. We thus decided to start from 
a more intrinsic completely random initial condition where N^/q up spins are randomly 
distributed on the lattice, and considered 800 x 800 square lattices, with equal nearest and 
next nearest neighbor couplings. For each value of g, the results have been averaged over 
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80 samples. These large sizes are necessary in order to allow large coarsening times. In 
fact, notice that the smaller is q, the shorter are the accessible coarsening times. Indeed, 
for a small concentration of minority phase (as for q = 100, 200), a too small number of 
droplets is obtained after a rather short time. 

We found that the scaling function g, once normalized such that g(0) = 1, is only very 
weakly q dependent, and is almost indistinguishable from the curves already presented on 
Fig. 4. More interestingly, we computed the autocorrelation A(t) = (sj(i)sj(O)) — uiq ~ 
L(£) _/ \ In order to determine A properly, and despite the long accessible coarsening 
times, we also had to use an interpolation scheme introduced in [9] for the usual Ising 
model. The authors argued that the effective exponent X(t) measured at time t should 
behave as, 



where b is a time scaling factor chosen in the range 10 — 40 depending on q and the 
speed of the simulation. Aoo is the exponent to be found. This relation was fairly 
well obtained for all q and typically modified the naive value of A (measured at large 
time) by 0.03 — 0.06 depending on q. We insist again on the fact that, due to the 
observed sensitivity to initial conditions (totally random, preformed bubbles...), it is 
possible that systematic errors are actually comparable or even bigger than the error 
bars [24]. The results of these simulations are presented on Tab. 1. Although for 
certain values of q the difference between the obtained values of A for the Potts model 
and the globally conserved Ising model is bigger than the error bars (in fact only for 
q = 50), the two exponents remain very close for all q. Actually, except for q = 200, we 
systematically have Xp tts > ^Ising by a typical amount of 0.02. We cannot conclude 
whether this tendency is real or simply results from systematic errors in the simulation 
and the various extrapolation schemes used. However, considering the already mentioned 
problems affecting the Ising simulation, but also the Potts simulation for large q, we 
cannot exclude that these exponents are strictly identical, as suggested by our theoretical 

23 




L(t) 



-1 



(VIA) 



analysis. Also notice that for q = 2, the exact global conservation of the magnetization 
does not affect the value of A, which is consistent with the results in [18] [19]. 

VII. Conclusion. 

In this paper, we have studied in detail the phase ordering process following a tem- 
perature quench in systems possessing, in general, q (q > 2) degenerate ordered phases at 
low temperature. We studied the dynamics of the g-state Potts model which accurately 
describes these systems. We have derived several exact results in different limits and 
obtained an important and interesting correspondence between the two-point correla- 
tion functions of the g-state Potts model and that of an Ising model evolving with fixed 
magnetization uiq = 2/q — 1. This analogy has been particularly useful in extending 
the Gaussian closure approximation developed for q = 2 to the case when q > 2, and 
the results obtained from this approximation agree very well with our direct numerical 
simulations. A note about GCA is worth mentioning at this point. It is well known 
[2] that up to now, the GCA fails to capture the essential dynamics in many situations 
such as the model B dynamics where the order parameter is locally conserved or even 
model A with long-range interactions [2]. However, for short range model A systems, 
such as nonconserved Ising model, GCA has given reasonably good answer especially for 
the autocorrelation exponent A. So, it is not surprising that for our system, which is also 
a short range model A system, it produces reasonably accurate values for the exponents 
\(q). 

As mentioned in the introduction, the correspondence to the Ising model has also 
an interesting experimental significance. We hope that this study will motivate future 
experimental work to measure quantities such as X(q), as it has already been done for 
q = 2. 
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FIGURE CAPTIONS 



Fig. 1: Late time configuration for q = 30 fields evolving according to Eq. (V.2). Note 
the bubble marked by an arrow, which is highly non isotropic, resulting from the 
coalescence of two bubbles with the same index. 

Fig. 2: Potts model in one dimension: (a) the three elementary dynamical flips of the central 
Potts spin with index a; (b) associating l's to the a-phase and O's to the other 
phases, one obtains an effective Glauber dynamics; (c) associating a particle A to 
each interface, the model is mapped on a reaction diffusion model as described in the 
text. 

Fig. 3: Distribution of bubble radii for d = 2 and model A, in the mean-field approximation 
(Eq. (III. 4)), compared to the Lifshitz-Slyozov result for the locally conserved order 
parameter case (model B) [16] [24]. Both distributions have been normalized and the 
r axis is scaled such that (r) = 1. 

Fig. 4: Comparison of the scaled equal-time correlation functions generated by (a) numerical 
simulation of the Potts model with q = oo (in fact q equals the initial number of 
bubbles Nq ~ 32000), (b) numerical integration of Eq. (V.2) with q — 50 (symbols 
have the size of the typical error bars), (c) the mean-field theory of section III (Eq. 
(III.7)), (d) the large q GCA analytical result of section V. 



Tab. 1: A for different q as computed from the Potts model simulation, for the globally 
conserved Ising model (see text), and by solving the eigenvalue equation Eq. (V.15) 
numerically. 1 see also [9]. 2 see also [15]. 
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Table 1 (Sire & Majumdar) 



bbbbacccc &=> bbbbbcccc or bbbbccccc 

bbbbabbbb e=> bbbbbbbbb 

bbbbaaaaa z=t> bbbbbaaaa or bbbbaaaaa 

000010000 « 000000000 

000010000 000000000 

000011111 « 000001111 or 000011111 

— AA— r=t> — -A— or — A— - 



— A— - r=t- — -A— or — A— - 
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